ABSTRACT
The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has triggered myriad efforts to understand the structure and dynamics of this complex pathogen. The spike glycoprotein of SARS-CoV-2 is a significant target for immunogens as it is the means by which the virus enters human cells, while simultaneously sporting mutations responsible for immune escape. These functional and escape processes are regulated by complex molecular-level interactions. Our study presents quantitative insights on domain and residue contributions to allosteric communication, immune evasion, and local- and global-level control of functions through the derivation of a weighted graph representation from all-atom MD simulations. Focusing on the ancestral form and the D614G-variant, we provide evidence of the utility of our approach by guiding the selection of a mutation that alters the spike's stability. Taken together, the network approach serves as a valuable tool to evaluate communication "hot-spots" in proteins to guide design of stable immunogens.
ABSTRACT
COVID-19 is a highly infectious respiratory disease caused by the novel coronavirus SARS-CoV-2. It has become a global pandemic and its frequent mutations may pose new challenges for vaccine design. During viral infection, the Spike RBD of SARS-CoV-2 binds the human host cell receptor ACE2, enabling the virus to enter the host cell. Both the Spike and ACE2 are densely glycosylated, and it is unclear how distinctive glycan types may modulate the interaction of RBD and ACE2. Detailed understanding of these determinants is key for the development of novel therapeutic strategies. To this end, we perform extensive all-atom simulations of the (i) RBD-ACE2 complex without glycans, (ii) RBD-ACE2 with oligomannose MAN9 glycans in ACE2, and (iii) RBD-ACE2 with complex FA2 glycans in ACE2. These simulations identify the key residues at the RBD-ACE2 interface that form contacts with higher probabilities, thus providing a quantitative evaluation that complements recent structural studies. Notably, we find that this RBD-ACE2 contact signature is not altered by the presence of different glycoforms, suggesting that RBD-ACE2 interaction is robust. Applying our simulated results, we illustrate how the recently prevalent N501Y mutation may alter specific interactions with host ACE2 that facilitate the virus-host binding. Furthermore, our simulations reveal how the glycan on Asn90 of ACE2 can play a distinct role in the binding and unbinding of RBD. Finally, an energetics analysis shows that MAN9 glycans on ACE2 decrease RBD-ACE2 affinity, while FA2 glycans lead to enhanced binding of the complex. Together, our results provide a more comprehensive picture of the detailed interplay between virus and human receptor, which is much needed for the discovery of effective treatments that aim at modulating the physical-chemical properties of this virus.
Subject(s)
Angiotensin-Converting Enzyme 2/chemistry , COVID-19/virology , Polysaccharides/chemistry , SARS-CoV-2/physiology , Spike Glycoprotein, Coronavirus/chemistry , Amino Acid Sequence , Binding Sites , Glycosylation , Host Microbial Interactions , Humans , Molecular Dynamics Simulation , Mutation , Protein Binding , Protein Conformation , Protein Interaction Domains and Motifs , Virus AttachmentABSTRACT
The COVID-19 (coronavirus disease 2019) pandemic underwent a rapid transition with the emergence of a dominant viral variant (from the "D-form" to the "G-form") that carried an amino acid substitution D614G in its "Spike" protein. The G-form is more infectious in vitro and is associated with increased viral loads in the upper airways. To gain insight into the molecular-level underpinnings of these characteristics, we used microsecond all-atom simulations. We show that changes in the protein energetics favor a higher population of infection-capable states in the G-form through release of asymmetry present in the D-form inter-protomer interactions. Thus, the increased infectivity of the G-form is likely due to a higher rate of profitable binding encounters with the host receptor. It is also predicted to be more neutralization sensitive owing to enhanced exposure of the receptor binding domain, a key target region for neutralizing antibodies. These results are critical for vaccine design.